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REVIEW OF INVESTIGATIONS PERFORMED IN THE U.S.S.R. ON CLOSE 
APPROACHES OF COMETS TO JUPITER AND THE EVOLUTION OF COMETARY 
ORBITS 

E. I. Kazimirchak-Polonskaya 


INTRODUCTION 

The problem of the cosmogony of the solar system has long occupied a 
central position in astronomical research. The origin of comets plays a 
substantial role in the solution to this problem. It is therefore necessary 
to study in detail the "original" and "future" orbits of the comets with 
nearly parabolic orbits and to examine in particular the great orbital 
changes that take place when comets pass within the spheres of action of the 
giant planets. 

In the U.S.S.R. — especially at the Institute for Theoretical Astronomy 
(I.T.A.) in Leningrad — considerable attention is paid to investigations of 
this type. In this respect it is useful to refer to lAU Symposium No. 45, 
organized in Leningrad in 1970, and in particular to the introductory report 
by Chebotarev (1972) . 

1. CALCULATION OF DEFINITIVE, ORIGINAL AND FUTURE 
ORBITS OF NEARLY PARABOLIC COMETS 

Early work in this field was carried out by Mikhajlov (1924), Sakk and 
Kulikov (1951), Dirikis (1953, 1954), Galibina (1953), and Shmakova (1953). 

Makover (1955a) developed a special method for calculating the original 
and future orbits of long-period comets; the method involves taking the true 
anomaly, rather than the time, as the independent variable. Several definitive, 
original and. future orbits have been determined (Dirikis 1956; Barteneva 1955, 
1965, 1970, 1971; Galibina 1953, 1963, 1964; Galibina and Barteneva 1965; 

Belous 1960,^ 1964, 1966, 1970). Galibina (1964) established that although the 
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overwhelming majority of the original orbits were elliptical, about one half 
of the comets with definitive hyperbolic orbits continue to have hyperbolic 
orbits in the future and will therefore leave the solar system- Very similar 
results were obtained by Brady (1965). Reference should also be made to 
Sekanina's (1966) general catalogue of definitive, original and future orbits. 

2. EARLY SOVIET INVESTIGATIONS OF THE GREAT TRANSFORMATIONS 
OF COMETARY ORBITS IN JUPITER'S SPHERE OF ACTION 

The first work in the U.S.S.R. involving the investigation of large 
perturbations on cometary orbits and the successful prediction of the returns 
of short-period comets was carried out by Dubyago and Lexin (1923) and by 
Dubyago (1924, 1925). 

Dubyago (1932a, 1932b, 1936, 1946, 1950, 1956a, 1956b) constructed a 
numerical theory for the motion of P/Brooks 2 from before the comet's discovery 
in 1889 to 1960, taking into account the perturbations by the planets Venus to 
Saturn and a variable secular acceleration. He also studied two passages of 
the comet through Jupiter's sphere of action, to minimum distances of 
Ajjiin = 0.000964 AU from Jupiter in 1886 and = 0.086 AU in 1922. In 

studying the first approach, he took into account the perturbations due to 
the Galilean satellites (although the effect proved to be negligible) and to 
Jupiter's oblateness. He also considered and dismissed the question, first 
raised by Poor (1894), of the possible collision of P/Broo'xs 2 with Jupiter's 
satellite V as the reason for the comet's disruption. 

Subsequently, Dubyago (1955a, 1955b, 1956c) was engaged in research on 
the motion of P/Shajn-Schaldach and the great transformation of its orbit 
that took place during a close approach to Jupiter shortly before the comet's 
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discovery in 1949. He also conducted theoretical investigations into the 
structure of comets and their possible disruption under the influence of 
Jupiter’s destructive forces (Dubyago 1942) and into the nongravitational 
forces that affect the motions of comets (Dubyago 1948). Some of his ideas 
on the structure of the cometary nucleus are closely related to those of IJhipple 
(1950, 1951); Dubyago (1948, 1956a) considered the dependence of the non- 
gravitational forces on solar activity, on the perihelion distance of the 
comet and on the orientation of the comet's orbit. 

An early work by the present author (Kazimirchak-Polonskaya 1950) included 
(1) a history of the studies of the motions of 32 short-period comets that 
approached Jupiter and other major planets; (2) a description of methods for 
considering the nongravitational effects on the motions of comets and a survey 
of the various hypotheses made between 1830 and 1950 on the causes of these 
effects; (3) the suggestion of a series of studies that might be made of orbital 
transformations for comets passing through and near Jupiter's sphere of action; 
and (4) the description of a new jovicentric method using special rectangular 
coordinates and taking into account the perturbations by the sun and planets, 
and the application of this method to a study of the motion of P/Wolf within 
Jupiter's sphere of action in 1922 (A^^n ~ 0.125 AU). The comparison of these 

calculations on P/Wolf with the observations in 1925, as well as with the 

/ 

calculations (using a heliocentric method of variation of arbitrary constants) 
by Kamienski and Bielicki (1935), was very favourable. 

Sochilina (1958) studied the changes in the orbit of P/du Toit-Neujmin- 
Dalporte when that comet passed near Jupiter's sphere of action in 1954 

“ 0.656 AU) and noted that the mean motions of the. comet and Jupiter 
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would then be very close to 2:1 commensurability; Fokin (1958) studied the 
extended passage of P/Oterma through Jupiter's sphere of action = 

0.168 AU) during 1936-1938, and Merzlyakova (1958) investigated that of 
P/Ashbrook- Jackson (A„j^j^ = 0.178 AU) in 1945. 

A re-examination by Kastel' (1965) of the very close approach of 
P/Brooks 2 to Jupiter in 1886 gave = 0.000985 AU, closely confirming 

the earlier result by Dubyago (1950). 

3. DIFFERENCES AND DIFFICULTIES IN METHODS USED FOR THE 
CALCULATION OF LARGE PERTURBATIONS BY JUPITER 

A statement of the problem and a review of research on the close approaches 
of short-period comets to Jupiter during 1770 - 1960 were given by Kazimirchak- 
Polonskaya (1961a, 1961b). The approaches of 33 short-period comets were 
discussed, and differences and difficulties in the methods used by the various 
authors were analyzed. The possibility of using these approaches to determine 
a more accurate value for the mass of Jupiter was demonstrated; such a 
determination has recently been made in the case of P/Wolf, for example 
(Kazimirchak- Polonskaya 1972a). A number of questions arise: 

(1) What method — variation of arbitrary constants, perturbations in 
rectangular coordinates, or whatever — provides the most accurate results in 
calculations on the transformations of cometary orbits in the spheres of action 
of Jupiter and other planets? A number of investigators (e.g., Rasmusen 1935; 
Herget 1947; Dubyago 1956a; Marsden 1963, 1967; Marsden and Schubart 1965; Stumpff 
1972; Klepczynski 1972) have applied Cowell's well-known method of perturbations 
in rectangular coordinates. On the other hand, Merton (1927), and especially 
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the Polish astronomers (Kamiehski 1925, 1926, 1948a, 1948b, 1951, 1957, 1959; 
BCamienski and Bielicki 1935, 1936; Kepinski 1958) have utilized the method 
of variation of elements, defending it as being the most accurate method for 
calculating passages through Jupiter's sphere of action. 

(2) . What kind of method — heliocentric or jovicentric — should be 
preferred for very deep penetrations of comets into Jupiter's sphere of 
action? This is a very cogent question, for there are often severe discrepancies 
between the results of heliocentric and jovicentric methods applied to the 

same calculation. 

(3) What differential formulae should be used in order to allow — without 
repeated integration — for small additional perturbations, such as those by 
Jupiter's satellites or by nongravitational forces? 

(4) What criterion should be used in choosing the step-size for the 
integration? 

4. THE CHOICE OF AN EXPERIMENTAL OBJECT FOR 
CHECKING THE VARIOUS METHODS 

In order to overcome some of the difficulties mentioned in the previous 
Section it is useful to select a special experimental object. P/Wolf is an 
appropriate choice for three reasons: (1) the numerical theory for the motion 

of this comet was very skillfully constructed by Kamienski (1959) for the two 
isolated intervals of time 1884 - 1918 and 1925 - 1959; (2) P/Wolf passed 
close enough to Jupiter in 1922 that the correctness of the calculations can 
be verified by examining the pre-1922 and post-1922 observations/ and (3) the 
nongravitational forces on P/Wolf are practically insignificant. 
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The planetocentric method has been developed in special coordinates 
(Kazimirchak-Polonskaya 1962a). Taking P/Wolf as an example, the author 
demonstrated the practical equivalence of the method of variation of arbitrary 
constants and the method in special rectangular coordinates, both in difference 
and in quadrature forms. Question (1) of the previous Section was therefore 
answered . 

Some advantages of the method in special coordinates, as opposed to 
Cowell's method, have been demonstrated (Kazimirchak-Polonskaya 1961c); the 
jovicentric form of the new method was worked out; and it has been demonstrated 
that the heliocentric and jovicentric methods give practically identical 
results in the case of P/Wolf. Table 74 of the cited paper contains the 
answer to Question (2) . 

Question (3) was solved by developing Encke's method in planetocentric 
form (Kazimirchak-Polonskaya 1962b) and by producing a series of differential 
formulae for taking into account various small perturbations. The procedure 
was applied to the calculation of the perturbations by Saturn on P/Wolf 
during the encounter with Jupiter in 1922. 

As regards Question (4), the author has developed a new criterion that 
gives the integration step size as a function of the distance of the comet 
from sun and all perturbing planets. 

5. CONSEQUENCES OF THE DEVELOPMENT OF ELECTRONIC COMPUTERS 

The rapid growth of electronic computers has opened up many new areas of 
research on the motions of the minor bodies of the, solar system (Kazimirchak- 
Polonskaya 1967a, 1967b, 1972b, 1972c; Kazimirchak-Polonskaya et al. 1958, 1972 
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Kazislirchak-Polonskaya and Terent'eva 1973). Among these numerous problems 
we ohall mention only two: (1) the construction of numerical theories of 

motion covering the whole period of obsen/ations of each comet, with full 
allowance for planetary perturbations and the effects of nongravitational 
forces; and (2) the investigation of the evolution of cometary orbits over 
the 400-year interval 1660 - 2060. The problems are closely related, and in 
practice the second one will be solved in conjunction with the first in the 
form of successive approximations. 

The remainder of this review will be concerned mainly with the second 
problem, which can be subdivided as follows; (a) studying the orbital 
evolution of short-period comets of two or more apparitions; (b) redetermining 
the orbits of the short-period comets of only one apparition and jthen investi- 
gating the orbital evolution by a special method; (c) classifying the various 
approaches to the major planets and establishing the principal characteristics 
of the evolution, of cometar3/ orbits; (d) studying the transformations of the 
orbits of fictitious comets passing within the spheres of action of Uranus and 
Neptune and examining the mechanism whereby comets may be captured by these 
planets; (e) elucidating of the role of the giant planets in the evolution of 
cometaxy orbits; (f) specifying the successive stages in the evolution of 
cometary orbits, with consideration given to the stellar perturbations and 
the diffusion theory for long-period comets; and (g) analyzing all hypotheses 
on the origin of comets and developing the most probable hypothesis linking, 
as far as possible, all comets into a single complex. 

OUIGMAL PAGE IS 
OF POOR QUALITY 
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6, INVESTIGATIONS ON THE EVOLUTION OF COMETARY ORBITS 
AND COMETARY CAPTURE 

In an important series of papers, Everhart (1967, 1968, 1969, 1970, 

1972a, 1972b, 1973) has applied and improved the statistical methods dating 
back to Newton (1878, 1893). Havnes (1972), using a simplified formulation, 
arrives at particular conclusions concerning the dominant influence of 
Jupiter on the evolution of cometary orbits. As a complement to these 
studies Kresak (1957, 1972a, 1972b, 1973) has made extensive investigations 
in which the Jacobi integral in the problem of three bodies is used for 
solving various cosmogonic questions, and somewhat similar approaches have 
been made by Vaghi (1973a, 19731)) and Lowrey (1973) • 

Marsden (1963, 1967, 1970) Marsden and Aksnes (1967), Stumpff (1972) and 
Klepczyncki (1972) have carried out exhaustive research on orbital transformations 
of different comets in Jupiter's sphere of action by Cowell's method, or more 
recently hy using the numerical integration program of Schubart and Stumpff 
(1966). 

The Soviet astronomers also invariably use numerical integration programs 
in theit research. Full allowance is made for planetary perturbations, and 
the methods are continually being improved in order to make them suitable for 
more and more precise modelling of real cometary motion, even when comets 
penetrate very deeply into the spheres of action of the major planets. 

7. METHODS AND COMPUTER PROGRAMS AT I.T.A. 

The complex of computer programs at I.T.A. includes routines for 
numerical integration, reduction of observations, comparison of calculations 


497 


with observations and improvement of orbits. At present there are in use 
three essentially different methods of integration and corresponding sets 
of programs for the BESM-4 computer; those by Belyaev (1972), Bokhan (1972) 
and Kazimirchak-Polonskaya (1967c, 1972b). They supplement and, if 
necessary, are used to check each other. 

In the set of programs by Belyaev the integration is performed by 
Cowell's method in single precision. Perturbations by Venus to Pluto are 
considered, and the step size ranges from 40 days to some tens of minutes. 

The choice of the step size is made automatically according to the criterion 
by Kulikov (1960). Nongravitational effects are not taken into account. 

In the Bokhan programs the method of variation of arbitrary constants 
by Herrick (1972) is used. Perturbations by MercuiTr to Pluto are included, 
and there is the possibility of allowing for nongravitational effects using 
the model by Makover (1955b) . The integration step is selected according to 
the criterion by Kazimirchak-Polonskaya (1967c, Table XII). The programs by 
Bokhan are intended mainly for investigating the motions of objects with highly 
eccentric orbits, notably P/Encke and (1566) Icarus. 

The programs by the present author are in double precision and take into 
account the perturbations by Mercury to Pluto and nongravitational effects. 

The choice of integration step ranges from 20 days to 5/64 day (1 hour 52.5 
minutes). The author’s set of programs, which includes some of the standard 
programs by Bokhan (1969, 1972), is especially suited for modelling the great 
transformations of cometary orbits in the sphere of action of any major planet. 
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8. THE USE OF METHODS AND SETS OF I.T.A. PROGRAMS 

FOR CONSTRUCTING NUMERICAL THEORIES OF COMETARY MOTION 

The Belyaev set of programs is the one that is in most widespread use, 
both at I.T.A. and at scientific centers in Kazan, Kiev, Tomsk and elsewhere. 
Among the comets investigated using these programs are P/Faye (Belyaev and 
Khanina 1972), P/Giacobini-Zinner (Evdokimov 1972), P/Tempel-Tuttle (Kondrat'eva 
1972), P/Stephan-Oterma (Shmakova 1972) and P/Ashbrook-Jackson (Merzlyakova 
1974) - The programs have also been applied to studies of the orbital stability 
of minor planets with "cometary" eccentricies (Chebotarev et al . 1970, 1972, 

1974) and of the orbital evolution of meteor streams. 

Using her own programs, the present author (Kazimirchak-Polonskaya 1972b) 
has eliminated the 1918 - 1925 discontinuity in Kamienski's (1959) theory of 
P/Wolf. Belous (1972, 1974a, 1974b) applied these same programs to P/Borrelly 
and to linking the two apparitions of P/Westphal and P/Brorsen-Metcalf (Belous 
1974c, 1974d). The two apparitions of P/Stephan-Oterma have been linked 
(Kazimirchak-Polonskaya and Belous 1974). Bokhan and Chemetenko (1974) have 
investigated the motion of P/Encke during 1901 - 1970, and Kazimirchak-Polonskaya 
and Terent'eva (1973) have investigated the motion and evolution of the orbits 
of various meteor streams. 

9. THE EVOLUTION OF COMETARY ORBITS DURING 1660-2060 

The orbital evolution during 1660-2060 has been studied for a total of 
52 short-period comets (although in a few cases, x^^hen an orbit was not 
sufficiently reliable, the interval was reduced to- only 200 years). The more 
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interesting results have been included in the following series of papers: 

Belyaev (1966, 1967, 1973a, 1973b); Belyaev and Khanina (1972); Belyaev and 
Reznikov (1973); Belyaev and Stal'bovskij (1974); Belyaev and Shaporev (1974); 
Belous (1974b); Kazimirchak-Polonskaya (1966, 1967a, 1967d, 1967e, 1967f, 

1971, 1972c, 1973). In the case of P/Wolf the effects of nongravitational 
forces were included too (Kazimirchak-Polonskaya 1967d), the comet’s variable 
secular deceleration having been changed after each approach to Jupiter 
according to the law established by Kamienski (1961) • Although there were 
two close approaches to Jupiter going back from the discovery date of 1884 
to 1660 = 0.12 and 0.25 AU) , comparison with a computer run in which 

the nongravitational forces were excluded shows that their influence was 
negligible. 

Belyaev. (1973a) studied the orbital evolution of P/Neujmin 2, first on 
the basis of the system of orbital elements by Neujmin (1948) , and then 
starting from elements he had determined himself. In spite of the 14 close 
approaches to Jupiter during the 400-year interval, the two results are very 
similar. 

Of course, one cannot generalize this finding to all comets: there are 

some rather exceptional comets (e.g., P/Lexell and P/Kearns-Kwee) where the 
smallest changes in the initial elements alter the whole course of the calculated 
orbital evolution quite dramatically. Other single-apparition comets whose 
orbital evolution has been studied are P/Schwassmann-Wachmann 3 (Belyaev and 
Shaporev 1974), P/Churyumov-Gerasimenko (Belyaev 1973b), and P/Gunn and 
P/Kojima (Belyaev and Reznikov 1973). 

^ POOR 
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10. THE CLASSIFICATION OF THE APPROACHES OF SHORT-PERIOD 
COMETS TO THE MAJOR PLANETS 

The study of the orbital evolution of 52 comets involved some 320 
close approaches to Jupiter, more than 40 approaches to Saturn and a few 
approaches to Uranus. There were 86 passages through Jupiter's sphere of 
action and one passage through that of Satuim. Close approaches of short- 
period comets to Jupiter and Saturn are certainly not infrequent events, 
and they follow a complex regularity the study of which is of interest from 
many points of view. 

On the basis of the available literature, Kazimirchak-Polonskaya (1967b) 
has classified 157 approaches of 63 short-period comets to Jupiter in terms 


11. THE PRINCIPAL CHARACTERISTICS OF THE EVOLUTION OF 

COMETARY ORBITS SUBJECTED TO GREAT PERTURBATIONS BY JUPITER 

The typical transformations of cometary orbits that arise as the result 
of passages through or near Jupiter’s sphere of action are illustrated in the 
following eight general examples: 

(1) . Comets that remain in Jupiter's family, but which are at first 
invisible from the earth, because their orbits have large perihelion distances 
and low eccentricities. After an approach to Jupiter, usually shortly before 
discovery, the perihelion distances are reduced, and the orbital eccentricities 
are increased. Examples are P/Churyumov-Gerasimenko (Fig. 1) and P/Gunn (Fig. 2). 

(2) . Comets that remain in Jupiter's family, but whose perihelion distances 
and orbital eccentricities pulsate in a generally irregular manner. 
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and the subscript 1-0 refers to Jupiter. The solid sections of the path of the comet 
are to the north of the ecliptic, the broken sections to the south. 
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Examples are P/Wolf (Fig. 3), P/Wolf -Harrington (Fig. 4) and P/Schwassmann- 
Wachmann 2 (Fig. 5). This type of motion was first described by Kamienski 
(1954). 

(3) . Comets captured by Jupiter from Saturn's family.* These comets 
initially had large perihelion distances, but successive approaches to 
Jupiter cause them to decrease until the comets could be discovered. In the 
course of time the perihelion distances may increase again, and the comets 
will become lost from view. Examples are P/Whipple (Fig. 6) and P/Comas Sola 
(Fig. 7). 

(4) . P/Oterma (Fig. 8) has an exceptionally unstable orbit. Jupiter 
captured it during 1936-1938 from Saturn's family into an orbit near 3:2 
commensurability. After an interval of some 20 years another approach to 
Jupiter caused it to be ejected back into Saturn's family. 

(5) . In contrast to P/Oterma^there must be^ comets that have stable orbits 
for extensive intervals of time. They are usually not observed and are 
located in the regions between neighbouring planetary families (Jupiter and 
Saturn, Saturn and Uranus, arid especially Uranus and Neptune), as well as 
beyond the orbit of Neptune. The one observable example of such a comet is 
P/Schwassmann-Wachmann 1, which is located entirely between the orbits of 
Jupiter and Saturn. 

(6) . Comets from the region between Saturn and Uranus initially having 
perihelia on the orbit of Jupiter and captured into the Jupiter family. 
Examples are P/Brooks 2 (Fig. 9) and P/Schwassmann-Wachmann 3, for which 

= 0.0056 AU in 1882 (Fig. 10). 


* We shall continue to define comet "families" on the basis of aphelion 
distance. However, it would be useful to consider at some time in the 
future the alternative definition proposed by Bielicki (1972). 
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Figure 4. P/Wolf-Harrington (Kazitnirchak-Polonskaya 1972c, p 
















(7) . P/Lexell (Fig. 11) was captured by Jupiter in 1767 0.018 AU) 

from a nearly circular orbit (perihelion distance q = 3.3 AU, period P = 10 years) 
into an elongated elliptical orbit with its perihelion inside the orbit of 
Venus (q = 0.67 AU, P = 5.6 years). The comet was discovered in 1770 as a 
bright object, owing to its exceptionally close approach to the earth (within 
0.016 AU) . The comet encountered Jupiter again in 1779, and this time the approach 
was so close (A ^y,-,- „ = 0.0015 AU) that the comet was ejected on an orbit having 

its aphelion far beyond the orbit of Pluto. The period of revolution increased 
to 260 years, and since the perihelion was removed to the orbit of Jupiter 
the comet will no longer be accessible to observation. As already noted, the 
future evolution of this comet is very sensitive to the initial conditions, 
and Fig. 12 shows the effect of changing the 1770 orbit slightly. It is not 
impossible that P/Lexell left the solar system on a strongly hyperbolic orbit. 

(8) . The orbital evolution of P/Keams-Kwee (Fig. 13) was investigated by 

us on the basis of the provisional elements determined by Marsden (1964) 
from observations covering an interval of six months. It illustrates fhe 
possible two-stage capture of the comet by Jupiter: during the first very 

close approach (A^^.^ = 0.042 AU) in 1855 Jupiter captured the comet from a 
hyperbolic orbit* into the Neptune family. After two revolutions around the 
sun, the comet passed deeply into Jupiter's sphere of action in 1961 (4^^^^.^ = 

0.032 AU) and was discovered two years later as a short-period comet of 
Jupiter's family. 


* It is necessary to stress that this is only a possible evolution and 
probably a mathematical fiction. A calculation by Marsden and Aksnes 
(1967), considering observations over an 18-month arc, indicates that 
the 1855 approach to Jupiter did not really occur. 
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12.. GREAT PERTURBATIONS OF COMETARY ORBITS BY SATURN 

We have also studied the orbital evolution of some of the short-period 
comets belonging to Saturn's family, notably P/Neujmin 1, P/Neujmin 3 and 
P/Gale. In the course of a 200-year interval P/Neujmin 3, for instance, 
experienced six approaches to Saturn, and four of them had quite considerable 
effects, causing perturbations of more than 10° in co and in 0. After a very 
close approach of this comet to Jupiter in 1850 = 0.12 AU) Saturn ceased 

to be the dominant influence on its evolution. During 1660 - 2060 P/Neujmin 1 
makes six approaches to Saturn and none to Jupiter, which shows that the 
secular evolution of the orbit of this comet is essentially determined by Saturn. 

Perhaps the most interesting comet of Saturn's family is P/Gale. 

According to the initial system of elements obtained by Dinwoodie (1959) from 
the two apparitions, between 1660 and 2060 P/Gale made nine approaches to 
Jupiter and eight to Saturn. One of the latter approaches is the only known 
passage of a comet through Saturn's sphere of action (\jxn “ 0.17 AU in 1798). 

As is happened, this approach did not result in a particularly great transfor- 
mation of the orbit of P/Gale. By varying Dinwoodie's elements slightly, we 
were able to decrease Ajjjx^ in 1798 to only 0.095 AU, and this caused a change of 
173° in 0. 

It thus appears that although Jupiter must play the dominant role in the 
evolution of cometary orbits, Saturn can exert a strong temporary influence 
on the evolution of the orbits of some of the comets belonging to its f ami ly. 


13. GREAT ORBITAL TRANSFOPvMATIONS OF FICTITIOUS COMETS 
IN THE SPHERES OF ACTION OF URANUS AND NEPTUNE 


Although we have investigated the orbital evolution of several of the 
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comets belonging to the families of Uranus and Neptune, (e.g., P/Stephan- 
Oterma, P /Pons-Brooks, P/Brorsen-Metcalf and P/Westphal) , there were only 
two minor approaches of P/Stephan-Oterma to Uranus, and none of the comets 
made any close approaches to Uranus or Neptune. In order to study the 
effects of such close encounters it was necessary to produce some fictitious 
comets. In order that these comets might bear some resemblance to real 
comets it is essential to discuss the existence of the Oort cometary cloud 
and the theory of diffusion and to consider the distribution of some of the 
orbital elements of long-period comets. 

Although the existence of the Oort cloud cannot be checked directly, 
some indication of its possible dimensions and structure may be determined 
by studying the "original" orbits of long-period comets for which fairly 
reliable definitive orbits have been calculated. The theory of diffusion, 
which dates back to van Woerkom (1948), has been elaborated by Oort (1950), 
Lyttleton (1953), Shtejns (1960, 1961, 1962, 1964, 1972), Shtejns and 
Kronkalne (1964, 1968), Shtejns and Riekst 3 rn'sh (1960), Shtejns and Sture 
(1962), Kendall (1961) and Whipple (1962). We know that comets that formerly 
belonged to the cometary cloud gradually diffuse into the inner part of the 
solar system. Consequently, there must be a concentration of invisible comets 
having low-eccentricity, low-inclination orbits and perihelia far from the 
earth. Many of them could be around the orbit of Neptune. 

The distribution of semimajor axes and perihelia of the orbits of long- 
period comets has been studied by many astronomers, notably Svedstrup (1883), 
Oppenheim (1924), Witkowski (1953, 1965, 1968, 1971, 1972) and Hurnlk (1959, 
1964) . They arrived at the conclusion that there .exists an incontestable 
connection between the distribution of the perihelia of these comets and the 
galactic equator, which fact testifies to the interstellar origin of comets. 
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Lyttleton (1948, 1953) suggested that the passage of the sun through a 
uniform interstellar dust cloud could lead to the formation of comets by 
accretion. 

In view of the above, we have considered in our studies comets having 
initial orbits of two types: 

(1) . Comets with nearly circular, low- inclination orbits that formerly 
belonged to the cometary cloud on the periphery of the solar system, but 
which have already approached the orbit of Neptune as the result of diffusion. 
These may also be the comets which belonged, according to a hypothesis by 
Whipple (1972) to an extensive belt of comets beyond Neptune’s orbit. 

(2) . Comets of cosmic origin arriving directly from interstellar space 

on hyperbolic orbits. These orbits are both direct and retrograde and penetrate 
Neptune's sphere of action. 

We shall designate the fictitious comets passing through the sphere of 
action of Uranus by U-1 and U-2 and those penetrating the sphere of action of 
Neptune by N-1, N-2, ..., N-8. In an earlier paper (Kazimirchak-Polonskaya 
1972c) we have treated the characteristic features of the orbital evolution 
of the fictitious comets U-1, U-2, N-1, N-2 and N-3. 

14. ORBITAL EVOLUTION OF COMETS N-4 TO N-8 
Comet N-5 has a nearly circular transplutonian orbit of small inclination 

O 

(i =8.7); its perihelion is located near the orbit of Neptune, and its 
revolution period is 210 years. Having penetrated deeply into Neptune's sphere 
of action in 1715, = 0.0004 AU) , this comet leaves along a direct orbit 

O 

(i = 22.1) having its perihelion between the orbits of Uranus and Saturn, its 
aphelion not far beyond the orbit of Neptune and a revolution period of 103 years 
(Fig. 14). 
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Figure lA. Comet N-5. 


Comet N-6 is of interstellar type; it has a retrograde hyperbolic 

O 

orbit (i = 135.2), and its perihelion is located between the orbits of 
Neptune and Uranus. Its encounter with Neptune occurs in 1710 = 0.0005 AU)j 

before perihelion passage. The comet is ejected along a nearly circular, direct, 
transplutonian orbit (i = 16?0) having its perihelion near the orbit of Neptune, 
its aphelion beyond the orbit of Pluto, and a revolution period of 220 years 
(Fig. 15). 

O 

Comet N-7 is another retrograde interstellar comet (i = 159.0), but its 
perihelion is located between the orbits of Uranus and Saturn. It encounters 
Neptune after perihelion, passes twice within a very small distance of the 
planet =? 0.00035 and 0.00065 AU) . Afterwards it retains its retrograde 

o 

orbit (i = 136.6), but with a perihelion distance of only 1.2 AU; its aphelion 
is located near Neptune's orbit, and the revolution period is 61 years (Fig. 16). 
The final orbital inclination of comet N-7 is very similar to that of P/Pons- 
Gainbart, while the size and shape of its orbit are practically identical with 
those of P/Westphal. 

Comet N-8 also has a hyperbolic retrograde orbit, but after penetrating 
Neptune's sphere of action, it leaves on an orbit that is even more hyperbolic 
than initially. 

Finally, we briefly mention Comet N-4, discussed in detail elsewhere 
(Kazimirchak-Polonskaya 1975). It is captured as a stable satellite of Neptune 
having an orbit intermediate between those of Triton and Nereid (Fig. 17). 

We conclude that the planets Jupiter to Neptune, with their great masses 
and extensive spheres of action, have a substantial effect on the evolution of 
the orbits of comets. These planets can transfer comets from one planetary 
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Figure 15. Comet N-6. 






family to another and in exceptional cases can remove them beyond the 
limits of the solar system; conversely, they can also capture comets from 
transplutonian and even from hyperbolic orbits. 


15. COMETARY ORIGIN 

In a number of papers (Kazimirchak-Polonskaya 1967a, 1967f, 1972c) the 
author has criticized the classical theory of cometary capture and demonstrated 
that the various simplifications made — though perhaps fully justified at the 
time — were the cause of contradictions that arose between theoretical 
deductions (Callandreau 1892; Newton 1893) and the observations. 

Using electronic computers and the modem methods of celestial mechanics, 
we reconsidered the numerical theory of capture and presented the successive 
stages of cometary orbital evolution, taking into account stellar perturbations 
and the theory of diffusion. This enabled us to eliminate the discrepancies 
that formerly existed. Our more recent investigations confirm our view that 
cometary capture undoubtedly takes place, but it is a very complex process, 
extending in some cases for perhaps millions of years. 

The process begins either with the sun capturing interstellar matter into 
its own very extended sphere of action, or with the formation, in some way, of 
a cometary cloud at the periphery of the solar system. The problems of the 
stability of cometary motion in the outer regions of the solar system and the 
shape of the cometary cloud have been studied by Chebotarev (1963, 1964, 1966, 
1970), Mezhinskij (1972), and Antonov and Latyshev (1972). 

Cometary perihelia are thrown into the inner part of the solar system, 
either as a result of the sun's capturing an interstellar comet or on account 
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of stellar perturbations on comets belonging to the cometary cloud. When 
such comets penetrate deeply into the sphere of action of one of the outer 
planets (particularly Jupiter or Neptune) they may be converted — immediately 
or in successive stages — into short-period comets. An interesting example 
of Jupiter's capture of a fictitious comet on a parabolic orbit has been 
discussed by Sitars'xi (1968) . 

The great majority of the comets in the cometary cloud will diffuse into 
the inner solar system in the course of 10^ or 10® years and, according to the 
diffusion laws, they acquire nearly circular orbits of small Inclination, with 
their perihelia concentrated at the distances of the outermost planets. 
Transplutonian orbits of small inclination and perihelia near the orbit of 
Neptune may be formed as the result of diffusion (e.g. , comet N-5) 
or by Neptune's capture of a comet from interstellar space (e.g., comet N-6). 

At this stage further captures can be made by the giant planets. These 
captures may develop by a slow evolution or catastrophically. Evolutionary 
capture consists of a series of successive steps, with capture first by 
Neptune, then by Uranus, by Saturn, and finally by Jupiter. Catastrophic 
capture may involve transformation from an interstellar orbit directly into 
the inner part of the solar system (e.g., comet N-7) or the capture by Jupiter 
of a comet directly from the families of Uranus or Neptune (e.g., P/Keams-Kwee) . 

This concept of capture permits us to suppose that all comets (whether 
interstellar or on nearly parabolic orbits, or of long or short period) 
represent a unified cometary system, the development of which is determined 
in the spheres of action of the sun and the planets Jupiter to Neptune. 
Nevertheless, it is desirable that there should be further critical analysis of 
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our numerical capture theory, as well as of other hypotheses of cometary 
origin (e.g., Lagrange 1812; Vsekhsvyatskij 1930, 1933, 1955, 1967, 1969, 
1972a, 1972b; Lyttleton 1953; Whipple 1972, Fesenkov 1972; Safronov 1972). 
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Table 1 


Short-period Comets of Neptune's Family 


Name and Designation of Comet 

j(A.U.) 

e 

Pfyr) 

i 

P/N-7 

1750 

1.23 

0.92 

61.6 

136° 

P ! Pons-Camlart 

1827 II 

0.81 

0.95 

63.8 

136 

P/ Westphal 

1852 IV 

1.25 

0.92 

61.2 

91 

Pj Duly ago 

1921 I 

1.12 

0.92 

67.0 
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Table 2 

Characteristic of Orbits of Neptune Satellites: Triton, Nereid and N-4 


Satel- 

lite 

Mean distance 
from Neptune 
( 10^ Km) 

^ min 

(lO^m) 


Sidereal period 
of revolution 

e 

of mean 
orbit 

i 

relative 
to ecliptic 

Triton 

359 

359 

359 

5^2&02’”V"' 

0.00 

132°. 79 

N-9 

529 

91 

957 

11 05 31 12 

0.73 

60.87 

Nereid 

5570 

1337 

9803 

35909 36 

0.76 

9.97 
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